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ABSTRACT 

The extraordinarily long outburst of GRS 1915+105 makes it one of the most remarkable low- 
mass X-ray binaries (LMXBs). It has been in a state of constant outburst since its discovery in 
1992, an eruption which has persisted ^100 times longer than those of more typical LXMBs. 
The long orbital period of GRS 1915+105 implies that it contains large and massive accre- 
tion disc which is able to fuel its extreme outburst. In this paper, we address the longevity of 
the outburst and quiescence phases of GRS 1915+105 using Smooth Particle Hydrodynamics 
(SPH) simulations of its accretion disc through many outburst cycles. Our model is set in the 
two-a framework and includes the effects of the thermo- viscous instability, tidal torques, irra- 
diation by central X-rays and wind mass loss. We explore the model parameter space and the 
examine the impact of the various ingredients. We predict that the outburst of GRS 1915+105 
should last a minimum of 20 years and possibly up to ~ 100 years if X-ray irradiation is 
very significant. The predicted recurrence times are of the order of 10 4 years, making the 
X-ray duty cycle a few 0.1%. Such a low duty cycle may mean that GRS 1915+105 is not 
an anomaly among the more standard LMXBs and that many similar, but quiescent, systems 
could be present in the Galaxy. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXB) are bright X-ray sources con- 
sisting of an accreting black hole/neutron star primary and a Roche- 
lobe filling, low-mass secondary star. Matter is transferred from the 
secondary star via the inner-Lagrange (LI) point and forms an ac- 
cretion disc around the primary. In some cases the accretion disc 
undergoes sporadic outbursts which are thought to be triggered by 
a thermal-viscous instability, resul ting in an inc reased mass ac- 
cretion rate onto the primary (e.g. lLasotall200lh and a rapid X- 
ray brightening. Many transient systems have been ob served in 
the lo cal neighbourhood (Milky Way, LMC, SMC, e.g. iLiu et all 
and the duration of the outbursts is typically of the or- 
der of months. The quiescence period between outbursts in these 
systems generally lasts one to a few years, making their duty cy- 
cle ~ 1%. Such a value for the duty cycle, based only on these 
short period well-studied systems, is the one generally assumed in 
the models of X-ray lu minosity functions of nearby galaxies (e.g. 
iBelczvnski et alj20 04"). However, nothing suggests that it holds for 
the systems that have not been observed going through a full cycle. 
This is the case, for example, of some long-duration transients for 
which the quiescence phase duration, hence duty cycle, can only 
be indirectly estimated from the cooling curve of the neutron star 



primarjQ (e.g. lWiinands et al.|[200ll ; ICackett et al.ll2008h . Further- 
more, such a value of the duty cycle could be valid for the shorter 
period systems only but not for their longer period counterparts. It 
is this latter question we address in this work by exploring the case 
of the long period system GRS 1915. 

GRS 1915+105 i s one of the most well studied LMXBs 
(see Fender & Belloni (2004) for a review). Recent observations 
dHarlaftis & Grein er 2004) suggest the system contains a 14.0 ± 
4.4 Mq black hole with a 0.8 ± 0.5 Mp sec ondary, in a binary of 
orbital period 33.5 days jGreiner et afl|200lh. The outburst which 
lead to the system's detection in 1992 JCastro-Tirado et al ] 19921) is 
still proceeding to this day, lasting 20 times longer than any other 
LMXB. This longevity can be explained in simple terms by the 
long orbital period, which means there is plenty of room in the sys- 
tem for a very large ac cretion disc to form around the black hole 
(-Rdiac ~ 2 x 10 12 cm. lTruss & Donej|2006l) . A large disc ensures 
there is a large reservoir of mass available to fuel the long out- 
burst. Analytic estim ates of the dura t ion of GRS 1915+105's out- 
burst were derived bv lTruss & Pond d2006t ). By their nature, these 
estimates were not able to account for the dynamical and viscous 
evolution of the accretion disc. In this paper we re-address the dura- 
tion of the outbursts of GRS 1915+105 and also consider its X-ray 
duty cycle using smooth particle hydrodynamics (SPH) simulations 



1 For example, Iwiinands et alj feOOll) estimated that the quiescence phase 
of the quasi-persistent transient KS 1731-260 could be several hundreds 
E-mail: pd75@astro.le.ac.uk years long. 
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of the accretion disc through many outburst cycles. In section[2]we 
present the physical effects included in the simulations and the as- 
sociated free parameters. Section [5] deals with the implementation 
of these ingredients in the SPH code. In section [4] we present our 
results and estimate the duty cycle of GRS 1915+105. In the last 
section we discuss the implications of this work. 



2 PHYSICAL INGREDIENTS 

2.1 The disc instability model (DIM) 

The disc instability model was first developed to explain the 
outbursts of cataclysmic variables (dwarf novae) and has been 
successfully extended to explain the X-ray out bursts of LMXBs 
jMineshige & Wheelerll 19891 ; iDubus et ai]|200lh. He re, we briefly 
address the DIM and refer the reader to lLasotal d200ll) for a review. 

Vertical (the direction perpendicular to the orbital plane of the 
binary) thermal equilibrium in the disc results in a relation between 
the surface density of the disc E and its temperature T at any given 
radius. The locus of equilibrium positions in the E — T plane, takes 
the shape of an S-curi<e. The negatively sloped middle branch is 
caused by the sudden change of opacity as hydrogen ionises at 
T > 6500 K. When conditions dictate that the equilibrium posi- 
tion is situated on this middle branch the disc is forced to follow 
a limit-cycle between hot and cold states, defined by the critical 
surface densities E max and E m j n . These hot and cold states are 
thought to be respectively associated with enhanced and suppressed 
angular momentum transport within the disc. Many solutions utilise 
the S hakura-Sunyaev viscosity prescription dShakura & Svunvaevl 
Il973l) . namely, 



ac s H 



(1) 



where H is the scale height, c s is the local sound speed in the disc, 
and a a free parameter. The hot state is assumed to be associated 
with high disc viscosity au, and the cold state with low disc viscos- 
ity a c . This two-alpha approach was introduced early in the devel- 
opment of the models to reproduce observed outbursts durations. 
This translates into E max = E max (a c ) and E min = E m i n (ah). If 
anywhere in the disc E(r) > E max , the annulus enters the hot, 
high-viscous state, which propagate to nearby annuli. The front 
propagating inward forces the disc in the hot viscous state on its 
way: this high viscosity implies a high accretion rate onto the cen- 
tral object, leading to the X-ray outburst. The disc returns to quies- 
cence (low viscosity, small accretion rate) once E(r) < E m i n . 

It has been found that S max and E m j n scale (almost) linearly 
with radius l lCannizzo et aljl98^) : 
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and 



E min = 8.25flib M Mf "-"ai- -' 



gem 



(2) 



(3) 



where Mi is the primary mass in solar masses, and Rio is the radius 
in units of 10 10 cm. Equations (O and ([3} are used in our numerical 
setup, but note that some slightly different prescriptions exist. 



2.2 Irradiation 

As discussed above wherever E(r) > E max , the disc is in the hot 
highly-viscous state. However, this is not the only grounds for the 
disc entering the hot state. X-ray irradiation by the central object 
can keep the disc ionised and in the hot state out to a radius Ri rr - 



For typical black hole and disc parameters, iKing & Ritteri Jl998h 
find 



Rirr ~ 2.7 X 



M^cm, (4) 



where Mis is the central accretion rate in units of 10 18 g s _1 , r\ 
the accretion efficiency (we keep r\ = 0.1 throughout this work) 
and e the illumination efficiency. The latter parametrises the geo- 
metric properties of the disc (incident angle of irradiation, albedo), 
whereas 77 gives the fraction of the accretion energy that is radiated 
away. The accretion onto the black hole can, at most, proceed at the 
Eddington rate, 



„V0 _ ^O^Edd 



= 3.1 x 10~ 7 M yr" 1 



(5) 



for GRS 1915+105 (Ro is the black hole radius). Hence, from 
eq. 10, one gets the largest extent of the disc that can be switched 
to the hot state via accretion-powered radiation is 



R Edd = 1.21 x 10 



1/2 



(6) 



We cap the irradiation radius at its Eddington value so that Ri r 
Rvdd whenever Mi > Af Edd . 



2.3 Wind Loss 

Mass loss due to the local mass transfer rate exceeding the Edding- 
ton limit is also included in the model. The local accretion rate at 
radius r is given by, 



M(r) = -27rr?j r (r)E(r) 



(7) 



where v r (r) is the radial velocity. The Eddington rate at the same 
radius is. 



M Edd (r) 



rip 



GMi 



(8) 



where L Edd is the Eddington luminosity. If M(r) > MEdd(r) a 
wind should carry away the energy in excess of the Eddington limit. 
To allow for some freedom in the wind efficiency, we parametrise 
the latter by triggering the wind whenever 



M(r) > AM Edd (r) 



(9) 



The parameter A allows us to depart from the strict Eddington limit 
(A = 1): the smaller A, the smaller the threshold at which a wind is 
emitted. 

The three mechanisms detailed above contain the free param- 
eters that are varied throughout this work in order to reveal their 
relative effects, namely: ah, a c , e and A. Note also that the expres- 
sions for E max and E m i n are estimates only and that we allow for a 
change in their normalisations which gives us two extra parameters. 
This is detailed in the numerical setup below. 



3 NUMERICAL IMPLEMENTATION 

We study the long term ev olution of GRS 1915+105 using the 
SPH code first develo ped by [Murray] {l996|) and later modified by 
iTruss & Wvnnl d2004l) . It includes the thermal viscous instability, 
disc irradiation, and wind loss as described above. Moreover, the 
system is followed in the full binary potential, hence it naturally 
includes any tidal effects that may arise. 
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3.1 Particle injection and rejection 



3.3 Scaling outburst and recurrence timescales 



The particles are injected from the LI p oint, and into the primary's 
potential, with the transfer rate given bv lRitteil < fl99^ : 



Mo 



7.3 x 10" 



M 2 

M 



1 day 



M Q yr 1 



(10) 



Using GRS 1915+105 system characteristics, this gives — M2 ~ 
2 x KT 8 M yr _1 . 

Particles are removed from the simulation when they are 
within 0.04a (a is the binary separation) of the black hole, if they 
return to the secondary's Roche lobe or if they are at a distance 
r > a with a velocity greater than the escape velocity. The first 
condition implies that the accretion rates we derive are not the ac- 
cretion rates onto the black hole (as 0.04a 3> i?schw where Aschw 
is the Schwarzschild radius of the black hole). Some material may 
still be blown away in a wind before it reaches the black hole sur- 
face but resolution issues prevent us from studying the most inner 
regions of the accretion disc. 



3.2 Triggering the disc instability 

Practically, the disc is divided into one hundred annuli into which 
the surface density is evaluated. Whenever E(r) > E max (r), the 
corresponding ring is switched into the hot state (a = aj,). Con- 
versely, the disc switches back to the cold state (a — a c ) wherever 
E < S m in. Using Eqs. 10 and ([3} with the parameters of GRS 
1915+105 and the typical values aj, = 0.1 and a c = 0.01 , one 
gets 



E 
and 



max — i 1 max 



An 



2.5 x 10 



;cm (11) 



55 and = 10 are 

lax and A' m i n are also to 



Emm = A m in (-J ~ 2.2 X 10* (^-J g Cm" 2 . (12) 

In practice, we cannot use the values of A' max / min shown above. 
Building the disc up to A ma x ~ 2.5 x 10 5 gcm~ 2 would take 
a prohibitive amount of time: A^ix — K 
typical values we use. The quantities A' 
be varied to explore how a less or more massive disc would behave. 

Along the same lines, aj, and a c are also increased from their 
canonical values given above, in order to speed up the outburst 
and recurrence times. Doing so allows the system to quickly reach 
steady-state and undergo several outburst events during a single 
run. The draw-back is that results need to be scaled in order to get 
actual the recurrence and outburst times. This is detailed in $3.3\ 
below. 

The last technical point is that of the disc tra nsition between 
the co ld and hot states. It is performed following iTruss & Wvnnl 
J2004I) . and is set up to occur on the thermal time scale of the sys- 
tem, tth ~ fi^ 1 ' s me keplerian frequency). Once a switch is 
triggered, a follows 



where 



a(i) = a + ± a tanh ( it 



± _ (ah ± a c ) 



(13) 



(14) 



As discussed in Q3.2\ computational time considerations have 
forced us to use values for E max / min and a h / c which are differ- 
ent from their generally inferred values. In this section we detail 
the procedure we follow to scale the raw outburst and recurrence 
durations to that of the physical system. In the following, we de- 
note by real and sph the quantities related to the physical and 
the simulated system respectively. 

During an outburst most of the mass inside a certain radius, 
-Rout, will be accreted. At the start of the outburst, the surface den- 
sity inside this radius will be near E max (see Fig(3] left panel). The 
outburst timescale is then given by 



(15) 



where m is the accretion rate in the disc for r < R ou t and m is the 
total mass in that region. The latter is given by, 

/•-Rout pRout 

m~ / 27rrEmax(>)dr - / 27rrE m in(j")dr. (16) 
Jo Jo 

Using eqns. $2^ and <[3j, one gets 

m _ [Emax(jRout) _ Emin ( jRout )] . (17) 

3.05 

Assuming A ou t is constant, eqns. dl lb . J 1 2b and l |17| > give, 

m oc (A 

max Amin) ■ (18) 

Given that (see lFrank etall j2002h ) 

771 OC OfeEmax OC OhAmax (19) 



then 



to 



(A n 



A n 



A, 



(20) 



This provides a scaling for the results from the SPH code. The ratio 
of the physical outburst time to the simulated one is therefore, 



i real / rxrcal rxrcal \ sph r>-sph 

''out V max -^mirw^h -"-max 

,sph ~ I'ft'sph _ ft'sph rca l ivreal ' 
'out l^BKB ^minl^h -"-max 



(21) 



The quiescence time is scaled in a similar fashion; the time spent 
in quiescence is roughly the time needed to replenish the mass lost 
in the outburst by the mass transfer rate from the secondary — 7712 

t ■ ~— (22) 

—m,2 

where m is given by eq. J18t . Hence, 

A 

max 

-A 

m 1 11 



^quic 



-m,2 



(23) 



which gives the ratio of the physical to the simulated quiescence 
times as, 



AT' 
n 



s p h 



A 



sph 

max 



A. 



sph 



^sph 



(24) 



quicsc 

The recurrence time of the system is simply obtained by tree = 
tout + iqulesc- As will be seen later on, in the case of GRS 
1915+105, t q ui CS c > iout so that t rcc ~ tq U i csc . When interpreting 
the results based on these scalings, consideration needs to be given 
to the assumptions made deriving them. The main assumption is 
that Aout is independent of the value of E m ax- It is not exactly the 
case from one simulation to another, but stays true within a few 
percent. 

Our canonical simulation has the following parameters 
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^max — 55, min 

e = 7 x 10" 4 and A 



4.79, 



0.1, 



*f» = 1, a?> h = 0.1, 7? 

1.0. The values of a are typically a factor 



of ~ 

ar 1 



10 larger than the generally assumed values a c 
« 0.1. 



0.01 and 



4 RESULTS 

Outputs from the simulations are presented hereafter. The typical 
evolution of the accretion is described in ^4. II We then explore the 
parameter space of our simulations and report the effects on the 
disc outburst and quiescence durations in i )4.2ll4.3l and l4.4l Param- 
eters and outburst/recurrence times of all the simulations presented 
below are compiled in table [TJ 



4.1 Typical behaviour 



As stated in g3JJ particles are injected from the LI point and pro- 
ceed to build up in a disc around the accreting black hole. We 
present here a canonical simulation of GRS 1915's disc which was 
run using the parameters recorded in the field base of tableQ] The 
particle mass and K m3x are chosen to ensure approximately 10 5 
particles in the quiescent disc. We wait until the system reaches 
steady-state over several outburst^ to begin our analysis. It is im- 
portant to emphasise that we do not assume any given density pro- 
file for the disc but let it build up and find its own steady-state . 
This is in contrast to the analytic treatment o flTruss & Done] feOOfj) 
where the surface density profile of the disc was prescribed by the 
authors. 

The steady-state of the disc is illustrated in Fig. [Tj where the 
central accretion rate (top-panel), irradiation radius (second panel), 
wind loss (third panel) and number of particles in the simulations 
(bottom panel) are plotted against time. In the time interval de- 
picted in here, the disc goes through four outbursts during which 
the central accretion rate increases by at least three orders of mag- 
nitude up to ~ 10~ 5 M o yr -1 . The luminosity of GRS 1915+105 
makes it slightly super-Eddington, and the inferred typical accre- 
tion rate onto the black hole is ~ 10~ 7 Mq yr -1 . The discrep- 
ancy between these two values can be explained by the fact that we 
do not resolve the inner disc. The mass transfer of 1O _5 M0 yr -1 
in our simulations represents the transfer rate through the disc at 
n n = 0.04a 3> R sc hv The amount of mass crossing our inner 
boundary is not the amount of mass that will eventually reach the 
black hole and we expect most of it to be blown away in a wind. 
In fact, there are good reasons to expect that the mass transfer rate 
through such a large disc as is present in GRS 1915+105 must be 
much larger than ~ 10~ 7 Mq yr -1 and we discuss that point fur- 
ther in i]5.2l The accretion in GRS 1915+105 proceeds at super- 
Eddington rates and this explains the plateau observed in 7?j rr as 
the latter is capped at the Eddington radius (see t\2.2\ . This is a 
general feature of all our simulations. 

Caution is needed when reading outburst and quiescence du- 
rations from this figure as the plotted time is a direct output from 
the code: outburst and recurrence times have to be scaled according 
to ij3 .3 1 Doing so, one gets t ou t = 24 yrs and t qu i csc — 10722 yrs 
for this particular simulation. 



2 In this context, steady-state must be understood as the state reached by 
the disc when the number of particles oscillates around a constant average 
value when going through outburst and quiescence phases (see Fig.fT] lower 
panel). 




13000 14000 15000 16000 17000 18000 
time (days) 



Figure 1. From top to bottom: time evolution of the central accretion rate, 
irradiation radius, wind loss and number of particles in the simulation. 



Figures [2] and [3] follow the evolution of the surface density 
in the disc for the base simulation through an outburst. In fig. [3] 
the surface density at a given radius has been averaged over 27r. In 
each case, the three snapshots (a), (b) and (c) are taken just before, 
during, and after the outburst respectively. 

In fig. [3^, the surface density in the inner disc closely follows 
the E max (r) prescription. Exceeding this limit at r ~ 0.2a triggers 
the switch to the hot viscous state and marks the beginning of the 
outburst phase (see i j3.2b . The switch in viscosity is illustrated in 
fig. |4] showing the radial variation of the SPH a parameter at the 
start of the outburst. The particles in the annulus that enter the hot 
state spread to increase the surface density in neighbouring annuli 
and have them switch to the high viscosity state as well (fig. [3J3). 
Because of these matter waves propagating inward and outward, a 
significant portion of the disc enters the hot state, to be eventually 
accreted onto the black hole (but for the part blown away in the 
wind). 

Typically 20-30% of the particles (i.e. of the mass of the disc) 
are accreted during an outburst (Fig. [TJ lower panel). Once most 
of the material has been accreted, the surface density in the inner 
disc decreases to eventually become smaller than E m i n , as seen in 
fig. [3}:. This is the trigger for the disc to switch back to the cold 
low-viscosity state and ends the outburst phase. During quiescence 
the disc regains the mass it lost in outburst and the cycle repeats. 

Figure [2] is another representation of the above, where one 
clearly sees the inner disc being emptied throughout the outburst. 
The 2D plots have also the benefit to show the existence non axi- 
symmetric features, e.g. large spiral arms developing in the outer 
disc. 
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#* 
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Table 1. Table detailing each setup and results. The parameters used in all the simulations are given along with the raw and scaled outburst/quiescence times. 
The parameters where chosen so that Mi > — M% and the scaling method (see £|3.3t was applicable (for discussion see Sj4.2i . The setup i n alpha c 4 
never allowed the steady-state to be reached. The last three simulations high irr [1-3] explore large values of e for direct comparison to lTruss & Pond 
|2006|) in S|5jJ 

but could only be followed during one outburst, hence the lack of information on the quiescence time (see also discussion in §44}- 



4.2 Varying the viscosity 

The a-parametrisa tion of the viscosity in accre tion discs has been 
widely used since IShakura & Svunvaevl ( fl973h . However, to date 
some confusion still exists around the actual value(s) that a should 
take to describe th e angular momentum transport in accretion discs 
dKing et alJkoOTl) . In this section, we allow for a c and ah to de- 
part from their canonical values and report their effect on the disc's 
limit cycle. The following results have been obtained from the sim- 
ulations alpha h [1-4] and alpha c [1-4] which full set 
of parameters can be found in tab. Q] For numerical reasons, we 
could not change -K"^f ax / min with ot c /h as ec i s - HI 1 an d 0} require. 
Indeed, changing K Bph , . on a large range of values results in 

' & & max/ mm & & 

the loss of the transient behaviour we are studying here. Neverthe- 
less, the scaling of the outburst and recurrence times accounts for 
this technical detail. Figure [5] shows the effect of varying a/, (solid 
line) and a c (dashed line) on tout, tquicsc and duty cycle. Note that 
the values of a displayed are the physical values and not the SPH 
ones (cth/c = a h/c/10)' ^ e timescales are scaled as described in 
? 13 . 3 1 if not stated otherwise. 

The effects of changing a c are straightforward: little effect 
if any on the outburst time (fig. [5] top panel, dashed line) as one 
would expect and an increasing recurrence time with decreasing 
a c (middle panel). Reducing a c implies an increase of E ma x ac- 
cording to eq. l[2j. Assuming that the outer radius of the outburst 



time=15903.1 days, n = 125557 
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Figure 4. Radial profile of a at the beginning of the outburst phase. The 
inner disc is switching to the hot state while the outer regions remain in at 
low viscosity. 
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Figure 2. Snapshots of the surface density of the disc before (left), during (middle) and after (right) an outburst. The dashed circle represents the inner radius 
of the simulation. 



time = 15843.1 days, n = 124463 




time = 15903.1 days, n = 125557 




time = 16850.8 days, n = 93343 




Figure 3. Evolution of the azimuthally averaged surface density before (left), during (middle) and after (right) the outburst of fig. [2] Red and blue lines 
represent the critical surface densities S max and S m i n . 



and the density profile of the disc are only slightly affected by this 
changjfjthen more mass should be lost in the outburst for a smaller 
q c . The transfer rate from the secondary being fixed, this translates 
into longer quiescence times. 

The outburst time decreases with increasing cth as accretion 
proceeds much faster for larger values of the viscosity (fig. [5] top 
panel, solid line). As a result, the mass rate at the inner boundary 
increases making the flow even more super-Eddington. This how- 
ever, does not translate into more mass being lost in the outburst. 
The number of particles (not shown here) actually lost during an 
outburst decreases as au increases. With high values of ah, the in- 
ner disc empties quickly and quenches the outburst before the out- 
ward matter wave could turn more of the disc into the hot state. The 
outburst time is roughly governed by the viscous time calculated at 



3 This is verified in £|4.3I The trigger point of the outburst changes with 
the critical density profile but the radius to which the hot state propagates is 
roughly unchanged. 



the outermost part of the outburst R ou t and reads IPringld d 198 lh . 
1 



^out 



a h Sl K {H/RY 



(25) 



170 | ah ^ 1 ( ( H / R )° ut 



0.1 / \ 0.03 



Ro 



3/2 



yrs . 



where Sl/f = (GMi/R 3 ) 1 ' 2 has been used. The simulations give 
at most H/R ~ 0.03 and i? ou t ~ 0.3a, which gives 



I an yrs. 



(26) 



This simple analytical estimate is a rough upper limit for the out- 
burst time and is plotted as the dotted blue line in fig.[5](top panel) 
where it can be seen to reasonably reproduce the behaviour of the 
numerical results. 

The recurrence time also marginally decreases with increas- 
ing ah- This is understood when remembering that the quiescence 
phase of the disc is roughly the time needed to replenish the mass 
lost in the outburst at a rate —Ah. As mentioned earlier, more 
mass is lost in the outburst for smaller values of a/, which ex- 
plains the trend reported in fig. [5] (middle panel, solid line). Note 
that the two triangle points have not been scaled using the method 
reported in £13.31 Indeed, the latter is only valid if Mi 3> —M2. 
This would always be the case if working with the real disc mass, 
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Figure 5. t ou t , t rcc = t ou t + iquicsc and duty cycle shown as a function 
of a/j (solid line) and a c (dashed line). The two triangles represent simu- 
lations that have not been scaled by the method described in £|3.3I See text 
for details. 



but because of the numerical setup (small if max and Jf m i n ), we 
only get Mi > —M2 for the lowest values of an. This means 
that, for these two simulations, the disc significantly replenishes 
while still in outburst. This indicates the disc is close to persistent 
outburst. However, it is a limitation of the simulations and limits 
the range over which parameters are allowed to vary. The duration 
measured from the evolution plots (like that of fig. Q} largely un- 
derestimates the actual quiescence times. To overcome that issue, 
the quiescence times of these two points have been explicitly cal- 
culated from t qu i, 

lost in the outburst and not from the raw value of t a 



Mout / — M2, where M out is the total mass 

output 
from the code. 

The lower panel in fig.[5]shows the variation of the duty cycle. 
Both etc and an affect it in opposite ways, similar in amplitude. 
Nonetheless, it spans on a very narrow range with an average value 
from 0.1% to 0.45%. 



4.3 Changing the critical density profiles 

The critical density profiles we use are given by eq. i ll It and 
d 1 2b . which are fitted on the vertical disc structure obtained by 
ICannizzo et al.l dl988l) . Different fits can be obtained depending 
on the disc model (irradiated or no t) used and other prescrip- 
tions can be found in the literature 1 Cannizzo & Wheeler] 1 19841 ; 



lHameurv et al.lll998l:lbubus et al.ll200ll) . Our setup does not allow 
for irradiation feedback on the disc structure or changing £ max / min 
accordingly. As also mentioned in the previous section, we do not 
have much freedom to change their normalisations i\" max / m i n and 
still retain the transient behaviour of the system. Within the range 
of values where the transient behaviour is present, fig. [6] displays 
the effects of ii"max (solid line) and K^. (dashed line) parameters 
on the outburst and recurrence times and on the duty cycle (simu- 
lations sig max [1-3] and sig min [1-4] in table fT). 

Increasing Jf^jJ implies a steeper slope for the critical surface 
density lK n . As the outburst proceeds and the surface density in 
the inner disc diminishes, the latter falls below quicker for 



Figure 6. tout , Wee and duty cycle as a function of K max (solid line) and 
A" m ; n (dashed line). 

higher values of K^J^ (also see fig. fj};). The outburst time is then 
smaller for higher values of K^ in (fig. [6] top panel, dashed line). 
The mass lost in the outburst is decreased in a similar fashion (see 
eq. l |18t ) so that the time needed to replenish the disc also decreases 
with increasing (middle panel, dashed line). 

The behaviour of the outburst and recurrence times with vary- 
ing ifmax 1S l ess straightforward. First of all, a steeper slope of 
Emax means a smaller trigger radius for the outburst. This is shown 
in fig. |7J where the snapshots are taken just before the outburst. In 
the upper panel, K^tx has a low value: the outburst is triggered in 
the outer disc, at r ~ 0.38a, and will propagate inward (outside- 
in outburst). For a high value of (lower panel), the surface 

density crosses its critical value at r ~ 0.12a and the front will 
propagate outward, up to i? ou t ~ 0.3a (inside-out outburst, similar 
to that of fig. [3). In an inside-out outburst, the outer radius of the 
outburst does not appear to depend on K^t*- tn ' s 1S tne case f° r 
the base and sig max [2-3] simulations. In the latter, the hot 
state always propagates up to R ou t ~ 0.3a so according to eq. i25i , 
the outburst time should remain unchanged. This is seen in fig. [6] 
(top panel, solid line). The only deviation from the analytical up- 
per estimate is for the lowest value of iCnax where the outburst 
propagates outside-in and is triggered at Rout ~ 0.38a. Used in 
eq. l !25t . it gives an outburst duration much larger than the simu- 
lated one. This indicates that this analytical estimate breaks down 
for outside-in outbursts. With increasing -A n iax> more mass is lost 
in the outburst, regardless of its type, and the quiescence time in- 
creases accordingly (middle panel, solid line). 

The duty cycle decreases with both if max an d -^mtn' Again, 
it spans on a narrow range around 0.25%. 



4.4 Effect of the irradiation and wind efficiency 

The last two free parameters of our model are the irradiation and 
wind efficiency, e and A respectively. The former defines the irradi- 
ation radius given by eq. © and the latter the threshold mass rate 
for the wind to be activqj as established by eq. (|9). The irradiation 



The higher A the larger the threshold, hence the less efficient the wind. 
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time = 16303.0 days, n = 100237 
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time = 8365.4 days, n = 134588 




Figure 7. Disc density profiles just prior outburst (black solid line). Top: 
simulation with K n f^ = 40 g cm -2 . Bottom: simulation with -R"max 
hi; g cm -2 . If sp . is the same in both cases. The triggering radius of the 
outburst is small for larger if^iax. 



efficiency depends on the disc structural and albedo of the gas (see 
Fran k et al. I Goo!)), two quantities our model does not include and 
we treat e as a free parameter. There is evide nce for slow winds 
being emitted fro m accretion discs in LMXBs dMiller et alj[2006l : 
lUeda et alj|200^) . However, our code does not possess the ingre- 
dients for the disc to consistently emit a wind, given its thermody- 
namical properties. Varying A allows us to artificially give more or 
less importance to the wind without resolving its origin. Figure [8] 
illustrates the impact of varying these two parameters on the tem- 
poral evolution of the disc (simulations irr [1-4] and wind 
[1-4] from table [D . 

The efficiency of the wind has no noticeable effect on the 
outburst time (fig. [8] top panel, dashed line). Indeed, R ou t is not 
changed significantly whatever A so that the viscous time in the 
hot state remains the same. The rate at which mass is lost in the 
wind is roughly equal to the central mass accretion rate but these 
high rates lasts only during a small fraction of the outburst dura- 
tion (fig.Q] t0 P an d third panel): the mass lost in the wind is much 
smaller than that accreted during an outburst. Therefore, the recur- 
rence time does not show any significant change, as long as A is 
high enough. For the smallest value of A however, the threshold for 
the wind is so low that a wind is emitted during quiescence. It then 
takes longer to replenish the disc, hence the slight increase of t rec 
when A = 0.8. The main effect of the wind during outburst is to 
reduce the central accretion rate (not shown here). 

Changing the irradiation efficiency leads to much more strik- 
ing behaviour (solid lines, crosses and diamonds). The system al- 
ways presents super-Eddington accretion so that the irradiation ra- 
dius is capped at Eddington and simply given by eq. The effect 
of the irradiation is twofold: i) it maintains the irradiated region in 



The disc structure itself depends on the irradiation heating. 
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Figure 8. t ou t , t rcc and duty cycle shown as a function of e (solid line) and 
A (dashed line). Crosses are for simulations where the particle injection is 
on during the outburst whereas it is turn off for the diamond simulations. 
See text for details. 



the hot state, even if S(r) < E m i n , ii) it can switch external re- 
gions where S(r) < E max into the hot state. Both effects should 
tend to increase the outburst time. For low values, up to e ~ 10 -3 , 
Rhr < Rout and only the first effect is at play (top panel, solid line, 
crosses). That regime only affects a small amount of mass and the 
outburst duration plateaus for e < 10 -3 and increases only when 
e ~ 10 -3 . The scale of the graph makes it difficult to visualise and 
we refer the reader to table[T]for the actual values. 

To explore the second regime, e > 10~ 3 is required. Because 
of the scaled-down disc we simulate, pushing e to higher values 
leads to Mi ~ —Ma and leads the disc to be in a continuous out- 
burst state. To overcome that difficulty, we turn off the injection 
of particles during the outburst phase (M2 = 0) in order to get 
the outburst to come to an end (this corresponds to the simulations 
high irr [ 1-3 ] in tableQ}. This is a valid approach since GRS 
1915+105 has Mi 2> —M%. In fig.[8]these simulations are repre- 
sented by the diamond symbol. In that regime, Rout — Rirr is 
much larger and a significant amount of mass enters the hot state 
in the outer disc, resulting in much longer outbursts. Remembering 



that .R^r c< e i/2 and i out ~ t viBC cc i? 3 ^, one gets t ou 
This is in reasonable agreement with the behaviour followed by the 
three diamond simulations, where t ou t oc e ' 7 . For the highest val- 
ues of e the entire disc is irradiated and 80% of its mass is lost in 
the outburst. 

The recurrence time follows directly from there. The longer 
the outburst time, the more mass there is to replenish, hence the 
longer the recurrence time. 



,3/4 



5 DISCUSSION 

5.1 The outburst duration 



iTruss & Done] d2006l) analytically estimated the outburst duration 
of GRS 1915+105. We briefly outline the framework of their cal- 
culation before comparing to our results. In their work, the outburst 
time was calculated following 
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Figure 9. Comparison between this work and lTmss & D one 1 2006) 



tout ~ 777 — r- ^ ' 

(Mdiso) 

The surface density in the disc was taken to follow £ max in the 
inner 10% of the disc before flattening off at larger radii. Hence, 

/■flout 

M out = / 27r J RE max (0.1i? disc )dR, (28) 

J0.1R diec 

where i? ut was simply the irradiation radius (eq. {4]>) . The average 
mass transfer rate through the disc read 

(Mdisc) = (Mi) + (M W ind) - Mi. (29) 

Remembering that A/2 <C (Mi) and using Mi = MEdd, the in- 
clusion of (A/ W i n d} as free parameter resulted in a maximum and a 
minimum length of an outburst for a given irradiation efficiency 
Their maximum outburst time is obtained when (M-wind) = and 
the minimum one when (M W i n d) — {Mi) = M^dd- 

In fig. [9l are plotted the outburst times predicted by 
iTruss & DoneThooel) and our work, against the irradiation effi- 
ciency. For values o f e < 0.7 x 1 ~ 3 , o ur outburst durations are 
larger than those of ITruss & Pond d2006t) . This is due to the fact 
that their model is solely based on the irradiation to switch the disc 
into the hot state. Therefore, for small values of e their predicted 
outburst times vanishes to zero. As discussed in ^4.41 in the low 
e regime of our simulations, the small irradiation radius has the 
mere effect of keeping a small amount of mass in the hot state a bit 
longer than without irradiation. When e > 10~ 3 and the effect of 
the irradiation becomes significant, the two approaches are in fair 
agreement. The duration of the current outburst of GRS 1915+105 
in conjunction with these models may provide some evidence of 
the irradiation efficiency in this system. 

The black hole LMXB 4U 1755-33 is very reminiscent of 
GRS 1915+105. It went through a ~ 25 year long outburst that 

6 This cannot be directly compare to our set up, where we showed that 
the efficiency of the wind had little influence on the outburst time (fig. [8] 
top panel, dashed line). It is difficult to make an analogy between our two 
approaches, but we can understand this apparent discrepancy in the follow- 
ing way: in our simulation, the measured Mi in rj n ^> i?Schw knows of 
the existence of a wind (i.e. the stronger the wind, the smaller Mi ) so that 
overall (M w ; n( j} ~ const and t ou t ~ const whatever the wind. 



ended around 1995, has a long orbital period of 4 . 4 dayfl and also 
shows jet-like sign atures l lAngelini & WhiteteOoHlPark et alj2005l: 
iKaaret et alj|200q) . To the authors' knowledge, the masses of the 
binary components have not been determined yet, which exclude 
any direct attempt to dynamically model this system. However, its 
outburst duration lies at the lower end of the range we find GRS 
1915+105 which is reasonable given that its orbital period, hence 
accretion disc, is smaller. 



5.2 Central accretion rate and the need for a wind 

The accretion rate measured at the inner boundary of the disc 
(ri n = 0.04a ~ 10 11 cm) in the simulations yield M n ~ 
10 -5 M© yr -1 . As mentioned in t]4.1| this is much larger than 
the Eddington rate, at which the black hole is observed to accrete, 
i.e. Mi ~ 10" 7 M yr" 1 at R S chw = 4 x 10 6 cm. 

So far, the outburst of GRS 1915+105 has lasted 17 years. 
Using eq. l |25l > with H/R = 0.03 and = 0.1 gives a minimum 
outer radius R ou t ~ 2 x 10 12 cm. Assuming the surface density 
profile of the disc follows E max (given by eq. Jilt), a rough upper 
estimate of the mass lost in the outburst is 

/"flout 

M out = / 27rrE max (r)dr ~ 2 x 10~ 4 M Q . (30) 
Jo 

If the outburst were to stop in the near future, this mass implies a 
transfer rate through the disc Mi n ~ 10~ 5 Mq yr -1 , in agree- 
ment with the transfer rate we measure at the inner boundary of the 
simulations. 

Conversely, if one takes the observationnaly inferred Edding- 
ton rate (Mi ~ 10 -7 Mq yr -1 ) as the transfer rate through the 
disc and use t out = 17 yrs, one gets M ou t ~ 10~ 6 Mq. Assum- 
ing the same density profile gives i? ou t ~ 2 x 10 11 cm, which 
is incompatible with a viscous/outburst time of 17 years. Hence, it 
would seem a wind or outflow must carry away the excess mass 
and we suggest that most (> 90%) of the mass flow through the 
disc must be lost i n this way. 

We note that lUeda et all J2009h recently concluded in the ex- 
istence of a thermally and/or radiation-driven disc wind in GRS 
1915+105 from the absorption lines in its spectra and inferred a 
launching radius ~ 10 J i?schw This corresponds to the innermost 
regions our simulations. However, they estimate the mass loss rate 
in the wind to be ~ 10 -7 Mq yr -1 , of the same order that the 
accretion rate onto the black hole. However, such mass loss rates 
are difficult to measure and we suggest that the loss rate could be 
significantly higher. 



6 CONCLUSION 

In this work, we have studied the long term evolution of the GRS 
1915+105 LMXB using global SPH simulations of its accretion 
disc, which include the two-alpha disc instability model, irradiation 
and wind. The parameter space has been scanned with the limita- 
tion of retaining transient behaviour in the simulations. We showed 
that all parameters (but for the wind efficiency) can play a signifi- 
cant role in the determination of the outburst and recurrence times. 

In agreement with both intuition and previous work, the out- 
burst duration is affected by ah (the smaller an, the longer the 



Long orbital period compared to usual LMXBs. 
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outburst) but not by a c . The opposite effect is noted for the recur- 
rence time (the smaller a c , the longer the recurrence time and a 
very marginal effect of ah-) 

Increasing the slope of the critical density profiles is also im- 
portant: i) the outburst is shorter if E m i n is increased but only 
marginally affected by E max , ii) an increase of E m i n and E max 
respectively decreases and increases the recurrence time of the sys- 
tem. 

When the irradiation efficiency is small (e < 10 -3 ) the vari- 
ations of tout and tquicsc are non negligible but span on a rela- 
tively narrow range: t ont £ [20 — 40] years, t qu i csc £ [7000 — 
20000] years and the duty cycle DC £ [0.1 - 0.5]%. However, 
as soon as the irradiation becomes non-negligible in determining 
the outer radius of the hot state region of the disc, the outburst 
duration increases considerably, potentially reaching 130 years for 
e = 5 x 10~ 3 . 

In lTruss & Pond J2006h . the authors relied on the irradiation 
to set the outburst duration. For large values of the irradiation effi- 
ciency, our results are in agreement with their analytical modelling. 
However, for small values of e, the thermal-viscous instability is at 
the origin of the outburst and we show that GRS 1915+105's out- 
burst can be expected to last at least ~ 20 ± 5 yr. As the outburst 
began in 1992 this raises the possibility that the outburst could end 
in the next decade. If so this would indicate that the X-ray irradi- 
ation efficiency is small, e < 1.5 x 10~ 3 . If however the outburst 
persists any longer, the conclusion that significant fractions of the 
outer disc are being irradiated is unavoidable. 

The extraordinary duration of the outburst of GRS 1915+105 
(and possibly of 4U 1755-33) is due to its long orbital period result- 
ing in an unusually large amount of gas available in its accretion 
disk. These extraordinarily long outbursts of long period, transient 
LXMBs should not be confused with the persistent outbursts of 
short period systems such as Sco-Xl or 4U 1957+1 1, in which sta- 
ble, persistent accretion can be easily establish ed by irradiation at 
moderate mass transfer rates ( lKingetalJl993) . 

Except if the irradiation efficiency is high, the duty cycle of 
GRS 1915+105 is ~ 0.5%, at best. Our canonical values gives 
DC = 0.2%: these values are slightly smaller than the typ- 
ical duty-cycles measured in shorter periods LMXBs (~ 1%), 
or that generally assumed for populatio n models and X-ray lu- 
minosity functions of ne arby galaxies dBelczvnski et al.l 120041 : 
llvanova & Kalogera] l2006h . This suggests than one needs to be 
careful when considering long period systems (for which very lit- 
tle (long-term evolution) data exists) and we will consider a more 
generic approach to these objects in a forthcoming paper. 
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